Simultaneous existence of two spin-wave modes in ultrathin Fe/GaAs(001) films 
studied by Brillouin Light Scattering: experiment and theory 
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A double-peaked structure was observed in the in-situ Brillouin Light Scattering (BLS) spectra 
of a 6 A thick epitaxial Fe/GaAs(001) film for values of an external magnetic field H, applied along 
the hard in plane direction, lower than a critical value He — 0.9 kOe. This experimental finding is 
theoretically interpreted in terms of a model which assumes a non-homogeneous magnetic ground 
state characterized by the presence of perperpendicular up/down stripe domains. For such a ground 
state, two spin-wave modes, namely an acoustic and an optic mode, can exist. Upon increasing the 
field the magnetization tilts in the film plane, and for H > He the ground state is homogeneous, thus 
allowing the existence of just a single spin-wave mode. The frequencies of the two spin-wave modes 
were calculated and successfully compared with the experimental data. The field dependence of the 
intensities of the corresponding two peaks that are present in the BLS spectra was also estimated, 
providing further support to the above-mentioned interpretation. 

PACS numbers: 75.70.-1 
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I. INTRODUCTION 



Ferromagnetic-semiconductor heterostructures, like ul- 
trathin Fe/GaAs(001) films, have received considerable 
attentioniiSiSj^ for their potential technological applica- 
tions in new magnetoelectronic devices A sharp and 
well ordered interface, without any dead magnetic layer, 
can be obtained, as bcc Fe grows epitaxially on GaAs 
thanks to the small lattice mismatch (1.4%). Together 
with the expected cubic anisotropy of bulk bcc Fe, a 
strong in-plane uniaxial anisotropy has been found in 
ultrathin Fe/GaAs(001) films, resulting from the atom- 
ing bonding at the interfaced The evolution of the latter 
anisotropy with film thickness has been quantitatively 
analyzed by some of us^ in a thorough in-situ Brillouin 
scattering study of the dynamical magnetic properties of 
such films. 

The same system Fe/GaAs(001) has now been found 
to display a very interesting phenomenon, for small iron 
thickness {tpe — 6 A): below a critical field He — 0.9 
kOe in-situ BLS spectra show a "double-peak" structure, 
therefore revealing the existence of two spin- wave modes 
for H < He- This feature is not completely new: it was 
already observed^ in Co/Au(lll) films, for tco > 6ML 
and H < He ^ ?> kOe. However, the novelty of our con- 
tribution is twofold, (i) The experimental observation 
of the double-peak structure has been done in two dif- 
ferent samples and both upon increasing and decreasing 
the magnetic field. This confirms that the phenomenon 
can be well reproduced and rules out the possibility it 
may be due to metastability effects, (ii) We develop a 
theory which explains the field dependence of the ob- 
served spin-wave frequencies as well as the intensities of 



the corresponding BLS peaks. 

The starting point of our theory is that the observed 
splitting of the spin- wave spectrum into two modes is not 
compatible with a homogeneous ground state. For this 
kind of system, the simplest and most natural explana- 
tion is to assume a perpendicularly magnetized up/down 
domain structure (for such low values of film thick- 
ness, in-plane magnetized domains are not energetically 
favourediSiiiii^) . Another possibility, a two sublattices 
spin arrangement, can be readily disregarded because the 
splitting is only observable below a critical value of the 
field and because it is hardly applicable to an epitaxial 
iron film. Our assumption is therefore an up/down do- 
main ground state. With increasing the magnetic field 
H, the magnetization of each domain gradually tilts in 
the plane and, for H greater than a critical value He, the 
ground state is homogeneously magnetized in plane. Two 
branches are found for H < He and a single (uniform) 
mode for H > He- 

At present an in-situ high resolution mapping of the 
magnetization is outside the reach of conventional micro- 
scopic techniques. In the absence of detailed information 
on the actual domain structure and for the sake of gener- 
ality, we are going to assume a stripe domains structure. 
Such a one-dimensional model has the advantage that 
the frequencies of the spin-wave excitations can be easily 
evaluated and their field dependence can be reproduced 
for different in-plane directions. 

The format of the paper is as follows. In Section II the 
experimental details concerning the sample preparation 
and the BLS technique are summarized. In Section III 
the frequencies of the spin-wave excitations with respect 
to a non-homogeneous ground state with up/down stripe 
domains are calculated using the Landau-Lifshitz equa- 
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FIG. 1: In-situ Brillouin spectra taken from a tpe = 6A thick 
epitaxial Fe/GaAs(001) film for different values of the ex- 
ternal magnetic field H, applied within the film plane along 
[— f 10], the hard in-plane direction. Two peaks, indicated by 
arrows, are simultaneously observed for H — 0.6 and 0.7kOe. 



tions of motion. We also estimate the field dependence 
of the spin-wave intensity of the two modes in the frame- 
work of a classical macroscopic model which relates the 
scattered intensity to the magnetization-dependent per- 
mittivity tensor. In Section IV the experimental results 
are presented and compared with the predictions of the 
theoretical model. Finally, in Section V the conclusions 
are drawn. Details about the calculation of the spin- wave 
frequencies can be found in Appendix A. 



II. EXPERIMENT 

Following the previous investigation of Fc/GaAs(100) 
films with different thickness,? in this work we focus our 
attention on 6 A thick Fe films, because for this par- 
ticular thickness BLS spectra exhibit the double-peak 
structure described in the following. The reproducibility 
of the results was checked studying two different sam- 
ples with the same nominal thickness. The two iron 



films were grown on a Si-doped GaAs(OOl) single crys- 
tal in an ultra-high-vacuum (UHV) chamber, specially 
designed to allow in-situ BLS measurementsiii^ (back- 
ground pressure 3 x 10^^" mbar) at GHOST laboratory, 
Perugia University. Details about sample preparation 
and structural characterization can be found in Ref. |^ 
About 200 mW of monochromatic p-polarized light, from 
a solid state laser (532 nm line), was focused onto the 
sample surface using a camera objective of numerical 
aperture 2 and focal length 50 mm. The back-scattered 
light was analyzed by a Sandercock type (3-|-3)-pass tan- 
dem Fabry-Pcrot interferometer .-iS" The external dc mag- 
netic field was applied parallel to the surface of the film 
and perpendicular to the plane of incidence of light i.e. 
in the so-called Damon-Eshbach geometry. BLS mea- 
surements of the spin-wave frequency were performed in 
situ at room temperature as a function of both the in- 
tensity and the in-plane direction of the applied magnetic 
field. Typical BLS spectra for such a film, taken with the 
magnetic field applied along the hard in-plane direction 
[-110], are shown in Fig. 1. The presence of a double-peak 
structure is evident for applied field values 0.4 < H < 0.9 
kOe. In contrast, a single peak was observed in the BLS 
spectra measured for magnetic field applied along both 
the [100] and the [110] directions (not shown). 

III. THEORY 

As discussed in the Introduction, the observed splitting 
of the spin- wave spectrum into two modes is not compat- 
ible with a homogeneous ground state. In the absence of 
experimental data about the actual spin configuration, 
we assume a simplified model with perpendicularly mag- 
netized up/down domains in the shape of parallel stripes 
of infinite length along the in-plane field direction and 
with vanishing thickness of the wall between opposite 
domains: see Fig. 2. While more comphcated up/down 
(or possibly canted) domain patterns cannot in principle 
be ruled out, the former one-dimensional model for the 
non-homogeneous ground state presents the advantage 
that the frequencies of the spin-wave excitations can be 
easily evaluated using the theory of domain mode ferro- 
magnetic resonance (DMFMR),^ii2ii242S and analytical 
results can be obtained for field applied in plane along 
high symmetry directions. Such a simple model turns out 
to be a useful tool to reproduce the spin- wave behavior; 
however, one should not expect it to be entirely realistic 
or able to account for other properties of the system, like 
the domain wall structure, the spatial dependence of the 
demagnetization factors, as well as the domain width and 
its field dependence. 

The free energy per unit volume of the system in Fig. 2 
is given by the sum of various contributions 

G = Gh + Gdip + G2± + G211 (1) 

where Gh is the Zeeman energy term due to the external 
field, Gdip is the term due to demagnetization, while G2i_ 
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FIG. 2: One-dimensional up/down stripe domain structure 
and coordinate system used to define the domain variables: z 
denotes the normal to the film plane. 



and G211 are the terms due to uniaxial out-of- plane and 
in-plane anisotropies, respectively. One has^i 

H AI 

Gh = [sin (^1 cos (j)i + sin 62 cos (f)2] 



Gd»p = |Af^2[cos6ii -f COS6I2] 



+ |m27V^,[cOS0i -COS02]' 

-I- — M^A'j,j,[sin Q\ sin — sin 6*2 sin 

(2) 



K2 



G2I. = |i[cOs2 0i+COs2 02] 



where Ms is the saturation magnetization; H is the ex- 
ternal field applied within the film plane (a;?;) along the x 
direction; K2^ > is a uniaxial out-of-plane anisotropy 
that favours the up/down spin alignment along z, the 
normal to the film plane; iV^z is the demagnetization 
factor associated to the up/down domain structure: since 
the stripes are assumed to be parallel to the field direction 
{x axis), one has = and Nyy = l — N^z- In general, 
the static demagnetization factor Nzz is a function of the 
domain aspect ratio L/t (where L is the domain width 
and t the film thickness) and of the rotational permeabil- 
ity fi 
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For magnetic field applied along the hard 
in-plane direction, the in-plane anisotropy contribution 
to the free energy is written as 



^2\\ = ^ [sin^ Oi sin^ 0i + sin^ 6*2 sin^ ^2] (3) 



while for field applied along the easy in-plane direction 



- ^ [sin^ 9i cos^ + sin^ 62 cos^ (^2] ■ (4) 



In the former case the uniaxial in-plane anisotropy K2\\ > 
favours the y direction (perpendicular to the field and 
to the stripes), and in the latter case it favours the x 
direction (parallel to the field and to the stripes). 

In terms of the free energy parameters, we define 
the out-of-plane anisotropy field H2± = 2K2±/Ms, the 
in-plane field H2\\ = 2K2\\/Ms, and the dipolar field 
Hdip = 47r Ms. It is customary to introduce the qual- 
ity factor Q — H2±/Hdip as the ratio between the out- 
of-plane anisotropy field H2±, favouring the perpendic- 
ular direction (2), and the dipolar field Hdip, favour- 
ing the film plane (xy). In the case under study, we 
have Q < 1. Finally, the saturation field is defined as 

Hsat = H2± — HdipNzz^^^ 

In the following, the equilibrium values of the polar and 
aziniuthal angles, obtained by minimizing the free energy 
Eq. will be denoted by the suffix "e". The frequen- 
cies of the spin-wave excitations are evaluatedil^ by the 
Landau-Lifshitz equations of motion (see Appendix A for 
details). Two modes, denoted by lo^ (acoustic mode) 
and uj~ (optic mode), are found for H < and a single 
(uniformSS) mode for H > H^. It is worth noticing that, 
despite their names, neither of the two modes is fully 
in-phase or fully out-of-phase: their peculiar character 
is disclosed by the analysis of the eigenvectors associ- 
ated to the two modes. In fact, assuming the external 
magnetic field to be parallel to the x direction (along 
which the domains arc infinitely long), one finds that 
the acoustic mode with frequency uj'^ is characterized by 
dynamic fiuctuations such that mi{t) + mf (i) = and 
m\{t) + m2{t) ^ 0. This corresponds to an out-of-phase 
precession of the dynamic moments Mi and M2 in the 
direction parallel to the domain wall and an in-phase pre- 
cession perpendicular to the domain wall. In contrast, 
for the optic mode the precession parallel to the domain 
wall is in-phase and the precession perpendicular to the 
domain wall is out-of-phase.^^ 

Hereafter we give the expressions of the spin-wave 
frequencies^^ when the field is applied in plane along the 
hard axis, see Eq. (O, or the easy axis, see Eq. Q. In 
both cases one has (/)ie — (f>2e = 0^ 

• Case (i): is along the hard in-plane axis. For 
H < He ~ Hsat, the minimum free energy is obtained for 
6*26 = TT — 6*12 where sin 6*16 = H/Hsat- The frequencies 
of the acoustic and optic modes are 



(-+/7)^ = 



Jj2 _ ff^fl _ HdipNyy 
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(5) 



respectively (7 is the gyromagnetic factor). For H > 
He = Hsat, one has — ^2e = 7r/2: the stripe domain 
structure is wiped out, and the film is homogeneously in- 
plane magnetized. The optic mode (i^^) disappears and 
the acoustic one (w^) becomes the saturated, in-plane, 
uniform^'^ mode with frequency 



(^/7)2 ^[H- H211] [H - {H2i_ - Hd^p)] 



(6) 
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• Case (ii): H is along the easy in-plane axis. For 
H < Hc = Hsat ~ ^2|| J the ground state has = tt — Oi^ 
where single = H/Hc and the frequencies of the acoustic 
and optic spin-wave excitations are 



(7) 



For H > He = Hsat — H2\\, one has 6'ie = 6*26 = 7r/2 and 
the frequency of the uniform mode is 

{lu/jY =[H + H2\\] [H - [H2i_ - Hd,p) + ff2l|] (8) 



IV. RESULTS AND DISCUSSION 

A. Spin-wave frequencies 

In Fig. 3 the measured spin-wave frequencies are plot- 
ted as a function of the intensity of the in-plane applied 
magnetic field H. When the field is parallel to the hard 
in-plane direction [-110] (Fig. 3a), two spin-wave modes 
are observed for 0.4 < H < 0.9 kOe. No differences in 
frequency, within experimental error, are found upon in- 
creasing or decreasing the field intensity, thus ruling out 
the possibility of metastability effects. In contrast, a sin- 
gle spin-wave mode is observed for H > 0.9 kOe. When 
the field is applied along either the intermediate [100] 
or the easy [110] in-plane directions (see Fig. 3b and 3c, 
respectively), a single spin- wave mode is experimentally 
observed at any investigated field intensity. The full-line 
curves in Fig. 3 are the theoretical spin- wave frequencies. 
The Hamiltonian parameters for the calculations were de- 
duced from previous experimental work on Fe/GaAs(001) 
films^ where a detailed fit of the BLS data was made for 
a whole set of samples with different values of the iron 
thickness, ranging between tpe = A A and 100 A. As 
tpe increases, the out-of-plane single-ion anisotropy K2± 
strongly decreases, while the dipolar field Hdip = 47rMs, 
favouring in-plane magnetization, increases; moreover, a 
biaxial in-plane anisotropy, favouring the [100] and [010] 
crystallographic axes, gradually develops. For the sam- 
ple with ti?e = 6 A, only the data at fields high enough 
for the magnetization to be homogeneous and in plane, 
were used to obtain the fit. The dipolar field was esti- 
mated to be Hdip — 17kOe,^^ the out-of-plane anisotropy 
field H2± = 13.9kOe and the in-plane anisotropy field 
H2\\ = 0.9kOei2& Using these parameters, the demagne- 
tization factor Nzz was self-consistently calculated^! for 
different values of the domain aspect ratio L/t. The vari- 
ation of the domain size with the applied field intensity, 
and thus the corrections to the frequency due to such 
variations, were neglected as being of a higher order ^ii^ 
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FIG. 3: Full-line curves: frequencies of the acoustic {i^'^), op- 
tic {uJ~) and uniform (uj) mode, as a function of the in-plane 
field H, calculated using H2\\ ~ 0.9kOe, H2± ~ 13.9kOe, 
Hdip = 17kOe, N^^ = 0.76 {Hsat = 0.98kOe), for three differ- 
ent cases: (a) field applied parallel to the stripes and along 
[-110], the hard in-plane direction (see Eqs. ||KJ and ©); (b) 
field applied parallel to the stripes and along [100], the in- 
termediate in-plane direction; (c) field applied parallel to the 
stripes and along [110], the easy in-plane direction (see Eqs. 
||7|l and Open circles: experimental data. 



The best overall agreement between theory and ex- 
periment was found assuming for the static demagneti- 
zation factor the value N^z = 0.76. Such an assump- 
tion, although it corresponds to a probably too low as- 
pect ratio^i is nevertheless able to justify the presence 
of domains in the system in spite of the fact that, for the 
considered iron thickness tpe = 6A, one has H2_\_ < Hdip 
(quality factor Q < 1). In fact, for the onset of up/down 
stripe domains, the condition H2± — HdipNzz > has to 
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be satisfied. As a further support to the domain hypoth- 
esis, it is worth observing that in epitaxial Co/Pt mul- 
tilayers a perpendicular (up/down) stripe domain struc- 
ture was indeed experimentally observed at remanence by 
magnetic force microscopy, while torque magnetometry 
measurements, providing Q < 1, had suggested a prefer- 
ence for in-plane orientation of the magnetization: 

Note that the two modes are well observable only for 
the case of field applied along the hard in-plane direction 
(Fig. 3a). Otherwise (see Fig. 3b, c) one has a consider- 
able shrinking of the coexistence region of the two modes 
and moreover metastability phenomena are likely to oc- 
cur since the energy of the stripe domain ground state is 
close to that of a homogeneous in-plane configuration. 



B. Light scattering intensities 

For the experimental backscattering geometry {H in 
plane parallel to the x axis and scattering plane perpen- 
dicular to H) the incident light has p-polarization (the 
optical wavevector and the optical incident electric field 
E/ have only y and z components) while the scattered 
light has s-polarization. Then its intensity is proportional 
to the square modulus of the x component of the polar- 
ization P induced in the film^^i^^i^^i^^i'^^ 



TO"=(i) - KE^j sin cos Oe 
2G44MsE'j cos 6le(sin2 9^ - cos^ 6^) 

my (t) \2G44MsEy sin Be - KEf 

TO^(t) \KEyshi^6e 

2G44MsE'j sin6le(sin2 9^ - cos^ 6*, 



(9) 



where K and G44 denote the first- and second-order 
(complex) magneto-optic coupling coefficients, respec- 
tively. For the film with up/down domains, we assume 
that TO"(i) = m"(i) -I- 1122 (t) {a = x,y,x) since the size 
of the laser spot is much greater than the lateral size of 
the domains. ■'■^ Taking into account that 4>ie = and 
sin 6ie — H/ He = h and expressing the dynamic fluctu- 
ations of the magnetization in terms of the fluctuations 
of the angle coordinates A9^{t), A(j)^{t) defined in Ap- 
pendix A, we obtain 

47rP^ = Ae-{t) \ - KE^ h{l-h^) 
~ 2G44MsEf (1 - h'^){2h^ - 1) 



A0+(i) 



2G44MsEy h? - KEf h 
- KE^ 



2G4AMsEf h^i2h^ 



1) 



(10) 



The field dependence of the intensities of the two modes 
can now be estimated taking into account that the eigen- 



vector associated with the acoustic mode is character- 
ized by A9^{t) = A(j)~{t) = while the optic mode has 
A9+{t) = A(t)+{t) = 0. 

• Acoustic mode with frequency co'^ . The intensity 
I^{H) of the light scattered by the acoustic mode has 
a maximum for H He since the fluctuations become 
very large. For zero field the intensity is zero, I^{0) = 0, 
since A9~{t) = and the coefhcients of the angle 
fluctuations are zero for h — 0. For H +00, the in- 
tensity vanishes, I~^{H) — > 0, since, upon increasing H 
above He, the uj^ mode evolves into the uniform mode 
uj and the fluctuations progressively decrease. This be- 
havior for I^{H) is similar to that of a perpendicularly 
magnetized uniform film^ 

• Optic mode with frequency uj~ . For H = He, the 
intensity of the light scattered by the optic mode is zero 
I- {He) = 0, since A9+{t) = A(j)+{t) = and the coef- 
ficient of the A9~{t) angle fluctuation is zero for h — 1. 
For zero field, the intensity (0) can be finite provided 
that the second order magneto-optic coupling coefficient 
is nonzero, G44 ^ 0. 

The field dependence of the intensity of both the acous- 
tic and the optic modes, as deduced from Eq. (^UJ, is 
qualitatively confirmed by the experimental spectra in 
Fig. 1. The intensity of the former mode exhibits a neat 
maximum for field values sligtly lower than He ~ 0.9 kOe 
and then it vanishes as the field is reduced below about 
0.5 kOe. The optic mode intensity, instead, shows a min- 
imum approaching He, in agreement with the theoretical 
predictions. 



V. CONCLUSIONS 

In conclusion, we have shown that a double-peaked 
structure is displayed by the Brillouin Light Scattering 
spectra of Fe/GaAs(001) films with tpe = SA when 
the field is applied in plane along the hard axis and is 
smaller than a critical value He — 0.9 kOe. The exis- 
tence of two peaks in the BLS spectrum should be the 
general feature of a film with a perpendicular domain 
structure (it is irrelevant whether the magnetization is 
canted or not). This feature disappears when H > He 
and the magnetization lies in the film plane. The rea- 
son why the unravelling of such a two-peaked structure 
in the Brillouin light scattering spectra of ultrathin mag- 
netic films is so rare might well be that many conditions 
have to be simultaneously satisfied. In fact, the optic 
mode, with frequency uj~ , has enough intensity in an ap- 
preciable range of fields only if G44 7^ and He are not 
too small. In contrast, the acoustic mode, with frequency 
Lu^ , has more chances to be observed since its intensity, 
though always vanishing in the H —^ limit, is expected 
to increase as H increases and to reach a maximum just 
at He- Another stringent requirement for the simulta- 
neous observation of two modes is that the competing 
out-of-plane anisotropy field H2± and easy-plane dipo- 
lar anisotropy field H^ip are of comparable magnitude 



6 



and that H2i_ — HdipNzz > 0, so that a perpendicularly 
magnetized up/down domain structure is energetically 
favoured for H ^ 0^ . This seems just to be the case of 
the Fe/GaAs films with tpe — 6A. In fact, for higher Fe 
thickness, one has Hdip ^ H2±, so that a homogeneous 
in-plane magnetized ground state is realized, while, upon 
reducing the Fe thickness, one would expect Hdip ^ ^2_l 
and a single spin-wave mode to be excited with respect to 
a homogeneous, perpendicularly magnetized metastable 
state. 

We hope that the results of this paper can stim- 
ulate other experimental groups to directly visualize 
the domain pattern, e.g. using magnetic microscopy 
techniques as a function of the external magnetic field 
intensity. ■^^ This should be done in situ, because the mag- 
netic anisotropy is strongly affected by the presence of a 
protective overlayer, so that formation of magnetic do- 
mains can be prevented^ 
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APPENDIX A: SPIN- WAVE FREQUENCIES 

The static equilibrium configuration of the system is 
obtained by minimizing the free energy G, Eq. with 



respect to the polar and azimuthal variables while the 
frequencies of the spin-wave excitations are evaluated^ 
by the Landau-Lifshitz equations of motion 



de 



'1,2 



7 dQ 



Pi, 2 



7 dQ 



dt 



sin^i^2 9(/)i_2 ' dt sin6'i^2 96'i^2 



(Al) 



where Q = 2G/Mg. The small oscillations of the system 
in response to an external perturbation are obtained by 
expanding in a Taylor series about its equilibrium value 
Qe up to the second order. Next, assuming for the set of 
variables a harmonic time dependence with frequency lo 
and introducing the normal coordinates = A^i ± 

A02 and A(/>^ = A0i ± A(/i2 (where A denotes a small 
variation), the equations of motion can be rewritten in 
matrix form as 



(A2) 



A+ 


— iz 





B+ 1 
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A(j)+ 





B- 


A- 


— iz 




Ae- 


B+ 





iz 


C" 




. A<p- 



where z = u sin 9\e and 



we denote the second 



dXidXj 

partial derivatives of the free energy with respect to the 
angular variables {Xi = 0i,X2 = <pi,X3 = 62, X 4 = 4>2)- 
The frequencies of the normal modes are obtained by 
imposing the condition for nontriviality of the solutions 
of Eq. HA2|I . i.e. the vanishing of the matrix determinant. 

For H > He the ground state is homogeneously in- 
plane magnetized, (0ie = 02e = 7''/2 and N^z — 1) and 
it results that A+ = A^ ^ Gn, B+ = B~ = 0, and 
C+ = = G22, so that a single, uniform mode^^ is 
obtained, with frequency (oj/j) = GiiG22- 
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